THE 3D DIMER AND ISING PROBLEMS REVISITED 

MARTIN LOEBL AND LENKA ZDEBOROVA 

Abstract. We express the finite 3D Dimer partition function as a linear combination of determinants of 
oriented adjacency matrices, and the finite 3D Ising partition sum as a linear combination of products over 
aperiodic closed walks. The methodology we use is embedding of cubic lattice on 2D surfaces of large genus. 



1. Introduction 

The motivation for this work has been a renewed interest of theoretical physicists in basic statistical 
physics models, the Ising and the Dimer models. Let us mention two such recent promising works. The first 
one is a new relation (duality) of topological strings and the Dimer model, see e.g. |UKV| . The second one 
is a relation of the theory of discrete Riemann surfaces [M] and Theta functions with criticality of the Ising 

problem [XMV1 EMS EM2- 

In this paper we present new formulas for the finite 3D Dimer and 3D Ising problems. The formulas are 
obtained in a combinatorial way. We hope that a profitable relation with the work mentioned above may be 
found. 

A model we consider here is the Ising version of the Edwards-Anderson model. It can be described as 
follows. A coupling constant Jij is assigned to each bond of a given lattice graph G; the coupling 

constant characterizes the interaction between the particles represented by sites i and j. A physical state of 
the system is an assignment of spin &{ £ {+1,-1} to each site i. The Hamiltonian (or energy function) is 
defined as H{a) = — ^2^jy eE Jij<Ji<Jj, The distribution of physical states over all possible energy levels is 

encapsulated in the partition function Z(j3) = e~^ H ^ from which all fundamental physical quantities 
may be derived. 

We may reformulate the Ising problem in graph theoretic terms as follows. A graph is a pair G = (V, E) 
where V is a set of vertices and E is now the set of edges (not the energy). A graph with some regularity 
properties may be called a lattice graph. We associate with each edge e of G a weight w(e) and for a subset 
of edges A C E, w(A) will denote the sum of the weights w(e) associated with the edges in A. 

An even subgraph of a graph G = (V, E) is a set of edges U C E such that each vertex of V is incident 
with an even number of edges from U. The generating function of even subgraphs TL{G, x) equals the sum 
of a w W over all even subgraphs U of G. A classic relation between the Ising partition function and the 
generating function of even subgraphs of the same graph states that 

Z{(3) =2" J] cosh (Mi) W(G,tanh(/3J y )). 

A subset of edges P C E is called a perfect matching or dimer arrangement if each vertex belongs to exactly 
one element of P. The dimer partition function on graph G may be viewed as a polynomial V{G, a) which 
equals the sum of a w ^ p ' over all perfect matchings P of G. This polynomial is also called the generating 
function of perfect matchings. 

The generating functions of even subgraphs and perfect matchings may be defined in a more general way 
as follows: associate a variable x e with each edge e of graph G, let x(A) = JleeA Xe an( ^ ^ e -&- ^ ne generating 
function of perfect matchings be the sum of x(P), P perfect matching of G. All results introduced in this 
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paper also hold in this more general setting; however the presentation using weights rather than variables is 
perhaps more natural. 



Convention for the cubic lattice. This paper studies properties of finite cubic lattices. Let us now fix 

some notation for them. 




(i,i,D 



Figure 1. Illustration of the linear order 
of vertices on a lattice Q(3, 3, 4). 



Let mi, m 2 , k be positive integers, where k is even and 
mi and m 2 odd. The cubic lattice Q — Q(mim 2 k) is 
the following graph: 

The vertices are: 

V(xyz), x=l,...,mi,y = l,.. . , m 2 , z = l,...,k. 

The edges are: 

• the vertical edges v xyz = {V xyz , V xy ^ z+ ^}, 
z=l,...,k-l, 

• the width edges w xyz = {V xyz , V x ( y+1)z }, 

y = 1,... ,7712-1, 

• the horizontal edges h xyz = {V xyz , V^+i)^}, 
x = 1,. ... mi — 1. 



Definition 1.1. We now fix as in Fig. 2] a linear order '<' on the set of the vertices of Q as 

V(ll), V(12), . . . , V(lm 2 ), V(2m 2 ), V(2(m 2 - 1)), . . . , V(mim 2 ) 
where V{xy) = V(xyl), . . . , V(xyk) and V xy denotes the reversal of V xy . 



v 
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• • • • • • • are m 2 x k 



Definition 1.2. Let us now define a set of pairs of edges R(Q). 
We let H be the set of all horizontal edges. For each x < mi we 
let W x consist of all width edges with their x-coordinate equal to 
x. We call sets H and W x layers of Q. It is convenient to depict 
the edges of a layer as the horizontal edges of a square grid, where 
lower left corner is the smallest in the fixed ordering of vertices. 
The first column of the layer H is ordered according to the fixed 
ordering of vertices. Size of layer H is mi X km 2 , sizes of all W x 



11 If e is an edge of such a layer L then let h(L,e) (v(L, e) respec- 

p. n . , , . tively) denote the horizontal (vertical respectively) coordinate of 

. . . r , , the first vertex of e m L. lor a layer L we let R(L) be the set of 

even) in a given layer L are marked . . u + i, + i,/r /\ ^ i./r N /r \/ 

all possible e' {h{L, e') is odd) to & 1 r P& f ( e ' ° f f ^ \ hat /l ( L ' 6 ) < e \f 

, . . . , v(L,e ) < V(L,e) + 1, and h(L,e) is even while n(L,e ) is odd, 

dg in p3ir with 6 — — - 

example on Fig. |21 Finally we let -R(Q) to be the union of R{L) 

over all layers of Q. 



2. Statement of the main results: The Dimer Problem 

Our result here is a simplification of the main result of |Llj . The reader is heartily encouraged to check 
the original statement for a comparison. Before we write the main theorems let us introduce some necessary 
notation. 
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Determinants. Q is a bipartite graph, which means that its vertices may be partitioned into two sets V\, V2 
such that if e is an edge of Q then | e (~l Vx j — |e n V2I = 1. We have \V\\ — | V2 1 = miiri2k/2. Let Z be the 
square (|Vi| x matrix defined by Zy = if {i, j} is not any edge of Q. If e = {ij} is an edge of Q 
(i £ Vi_, j £ V2) then .Zy = if i < j (in fixed linear order) and Zij = -a""' if z > j. 

We will consider matrix Z with its rows and columns ordered in agreement with the fixed ordering and 
we will assume that Vm £ V\. A signing of a matrix is obtained by multiplying some of the entries of the 
matrix by —1. 

An orientation of a graph G = (V, E) is a digraph D = (V, A) obtained from G by assigning an orientation 
to each edge of G, i.e. by ordering the elements of each edge of G. The elements of A are called arcs. We 
say that signing Z of Z corresponds to orientation D of Q if = —Z^ for (ij) £ 13 and i > j. For arc (ij) 
of orientation D of Q we let sign(D. (ij)) = 1 if i < j, and sign(D, (ij)) = —1 otherwise. 

Theorem 1. The dimer partition function V (Q . a) is 

V(Q,a) = 2~ c ^(-l)l« e ' e ') ei W)- si s n ( D ' e ) =si s n ( jD ^') = - 1 >ldet(Z(L>)), 

D 

where the sum is over all orientations of Q with all vertical edges positive, and C = km\(mi — l)/2 + 
kmi(m\ — l)/2. 

Knowing the proof Theorem ^ can be easily rewritten as 
Theorem 2. 

V(Q.a) = 2 c a w(M) - (2 C + l)£, 

where M is the unique perfect matching of Q consisting of vertical edges only and £ equals the average of 
dct(Z(D)) over all orientations D of Q satisfying: |{(e, e ) £ R(Q)', sign(Z3,e) = sign(D,e') = — 1}| is even 
(and again all vertical edges are positive). 

There are 2 2C orientations to be sum over in Theorem ^ and 2 C_1 (2 C + 1) terms in Theorem |21 to be 
averaged over. Therefore the exact numerical analysis is probably not possible. We just tested Theorem ^ 
for small planar lattices. However it could be interesting to study statistical properties of determinants 
det(Z(D)). Let us remark here also that Theorem ^ for a planar graph (2D Ising model, mi = 1) is still 
nontrivial, a question remains how to reduce it to the Kasteleyn's one determinant. And would it be possible 
to apply similar reduction even for the non-planar graph? 

3. Statement of the main results: The Ising Problem 

The theory described in section 0] combined with the results of |51 lL2| yields an expression of the 3D Ising 
partition function as a linear combination of products over aperiodic closed walks on Q. We get a remarkably 
simple formula, analogous to the statement for the dimer partition function. 

Products over aperiodic closed walks. Let G = (V, E) be a planar graph embedded in the plane and 
for each edge e let x e be an associate variable. Let A = (V,A(G)) be an arbitrary orientation of G. If 
e G E then a e will denote the orientation of e in A(G) and a" 1 will be the reversed directed edge to a. 
We let x ac = x -1 = x e . A circular sequence p — v%, a\, 1)2, d2> •••) (Vn+i — v±) is called non-periodic 
closed walk if the following conditions are satisfied: £ {a e , a^ 1 : e £ E}, <n 7^ a"^ and (a%, a n ) ^ Z m 
for some sequence Z and m > 1. We let X(p) — Y\"=iX ai . We further let sign(p) = (— 1) 1 +™(p), where 
n(p) is a rotation number of p, i.e. the number of integral revolutions of the tangent vector. Finally let 
W(p)=sign(p)X(p). 

Now assume p is a closed aperiodic walk in 3D cubic lattice Q. Can we define its rotation number? The 
following construction provides a solution. 

Assume some vertices of planar graph G embedded in the plane of degree 4 are marked. Then we call a 
walk p correct if it satisfies the 'crossover condition' at each marked vertex, i.e. it never enters and exits a 
marked vertex in a pair of neighbouring edges along the vertex. 

We naturally get a planar graph from Q if we consider Q drawn as bold Q(3,3,2) in Figure |3] and put 
a new vertex to each edge-crossing. These new vertices will form the set of the marked vertices. Clearly 
there is a natural bijection which associates to each walk p of Q the corresponding correct walk p' in the new 
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graph and we define rotation of p to be equal to rotation of p 1 . This defines W(p) for each closed aperiodic 
walk p in Q. 

There is a natural equivalence on non-periodic closed walks: p is equivalent with reversed p. Each 
equivalence class has two elements and will be denoted by [p]. We let H^([p]) = W(p) and note that this 
definition is correct since equivalent walks have the same sign. 

We denote by Y\(l + the formal infinite product of (1 + W([p]) over all equivalence classes of 

non-periodic closed walks of G. 

Let D be an orientation of Q. We further let 

na + w(\p]) = + w([ P ])\ Xc=sign(D ^ . 

D 

Theorem 3. Generating function of even subgraphs is 

H(Q, a) = 2-° J2(-l) l{{e ' e ' KR{Q) '' si S n (A e )=sign(.EX e ') = -l}| JJ (1 + W(\p\), 
D D 

where the sum is over the orientations of Q with all the vertical edges positive. 

Remark. Theorem [3] expresses the 3D Ising partition function as a linear combination of infinite products 
over aperiodic closed walks. For planar lattice this reduces to just one infinite product (conjectured by 
Feynman, proved in [5]) which provides an equivalence of planar Ising model with the quantum field theory 
of free fcrmion. There have been several attempts to generalize this equivalence for 3D. Earlier papers (e.g. 
[D*|1 attempt to replace closed aperiodic walks by 2D surfaces. Recent promising development mentioned in 
the introduction deals with discrete Remann surfaces and the theta function: 

It is known that the planar Ising partition function in the thermodynamical limit behaves like a 'common 
term' times a Riemann theta function. For torus, |AMVj asserts that in the thermodynamic limit and near 
to criticality, the Ising partition function behaves like a 'common term' times a linear combination of four 
Remann theta functions corresponding to torus. In Jm| and |CM1I ICM2 . there is evidence that for the 
critical Ising model, the dependence of the determinants of adjacency matrices on the Kasteleyn orientations 
is exactly the same as the dependence of the determinants of the Dirac operator, of the corresponding 
conformal field theory, on the spin structures of the 2D Riemann surface; this is given in terms of theta 
functions of half-integer characteristics. 

We believe that theorem [3] may provide a new insight into these efforts. 

4. Drawing of cubic lattices 
In this section we describe how we draw cubic lattices on 2-dimensional surfaces. 
4.1. Theory of generalised g-graphs. 

Definition 4.1. A surface polygonal representation S g of an orientable 2D surface of genus g consists of a 
base Bq and 2g bridges B^, i — 1, g and j = 1, 2, where 

i) Bq is a convex 4g-gon with vertices a±, a 4s numbered clockwise; 

ii) B\, i = 1, . . . ,<?, is a 4-gon with vertices x\ , x\ , x l 3 , x\ numbered clockwise. It is glued with Bq so 
that the edge [x\, x l 2 ] of B\ is identified with the edge [a 4 ( i _ 1 ) +1 , a 4 ( i _ 1 ) +2 ] of Bq and the edge [x\, x\] 
of B\ is identified with the edge [a4(i_i)+3, a 4 ( i _ 1 ) +4 ] of Bq; 

iii) B\, i = 1, . . . , g, is a 4-gon with vertices y\,y\, 2/3, y\ numbered clockwise. It is glued with Bq so that 
the edge [2/1,2/2] °f ^2 i s identified with the edge [a 4 (i_i) +2 , a 4 ( i _ 1 ) +3 ] of Bq and the edge [j/3, 2/4] of 
B\ is identified with the edge [a 4 (i_i)+ 4 , (Wi-ij+s] of Bq. Indexing is mod4g. 

Definition 4.2. A graph G is called a g-graph if it is embedded on S g so that all the vertices belong to the 
base B , and the embedding of each edge uses at most one bridge. We denote the set of the edges embedded 
entirely on the base by E and the set of the edges embedded on each bridge Bj by Ej, i = 1, . . . , g, j = 1,2. 
Moreover the following conditions need to be satisfied. 

1. the outer face of Go = (V, Eq) is a cycle, and it is embedded on the boundary of Bq, 

2. if e G E\ then e is embedded entirely on B\ with one end- vertex belonging to [xi,^] an( i the other 
one to [0:3, a: 4 ]. Analogously for e G E\. 
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We need a generalization of the notion of a <?-graph. 

Definition 4.3. Any graph G obtained by the following construction will be called generalized g-graph. 

1. Let g = gi + ... + g n be a partition of g into positive integers. 

2. Let S gi be a polygonal representation of a surface of genus gi, i = 1, ...,n. Let us denote the basis 
and the bridges of S gi by B l and B 1 - k , i = 1, n, j = 1, ...,<?; and fc = 1,2. 

3. For i = 1, n let be a <?;-graph with the property that the subgraph of iJ; embedded on Bq is a 
cycle, embedded on the boundary of Bq. Let us denote it by G\ 

4. Let Go be a 2-connected graph properly embedded on the plane and let F\, F n be a subset of the 
faces of Go- Let K l be the cycle bounding Fi, i = 1, ...,n. Let each K % be isomorphic to C l . 

5. Then G is obtained by gluing the H^s into Go so that each K l is identified with Ci. 

Orientations. Let G be a g-graph and let G* = (V, Eq U E % - ). An orientation Dq of Go such that each inner 
face of each 2-connected component of Go is clockwise odd in Dq is called a basic orientation of Go- Note 
that a basic orientation always exists for a planar graph. Further we define the orientation Dj of each G* as 
follows: Wc consider G* embedded on the plane by the planar projection of Ej outside B , and complete the 
basic orientation Do of Go to an orientation of G* so that each inner face of each 2-connected component of 
Gj is clockwise odd. The orientation — Dj is defined by complete reversing of the orientation Dj of G*. 
Observe that after fixing a basic orientation Dq, the orientation Dj is uniquely determined for each 

Definition 4.4. Let G be a g-graph, g > 1. An orientation D of G which equals the basic orientation D on 
Go and which equals D* or — on Ej is called relevant. We define its type r(D) G {+1, — l} 29 as follows: 
For i = 0, . . . ,g — 1 and j = 1, 2, r(D) 2 i+j equals +1 or —1 according to the sign of Dj +1 in D. 

Moreover we let c(r{D)) equal the product of Cj, % = 0, ...,g — 1, where q = c(r2j+i, ^21+2) and c(l, 1) = 
c(l, -1) = c(-l, 1) = 1/2 and c(-l, -1) = -1/2. 

Observe that c(r(D)) = (— 1)"2~ 9 , where n = |{i;r 2 i+i = ^+2 = — 1}|. 

For each generalized g-graph G we can define 4 9 relevant orientations D\, D^g with respect to a fixed 
basic orientation of Go, and coefficients c(r(Di)), i = 1, n in the same way as for a g-graph. Now we write 
theorem proved in |OL| which will be essential in proof of Theorem ^ 

Theorem 4. Let G be a generalized g-graph with a perfect matching Mq of Go- Let Dq be a basic orientation 
of Go. If we order the vertexes of G so that s(Dq,Mq) is positive then 

V(G,a) = Y,c(r(Di))PfG(Di,a), 

i=l 

where Di are all the relevant orientations of G. 

Definition of s(Do, Mo) follows from definition of Pfafjian. 

Let G = (V, E) be a graph with 2n vertices and D an orientation of G. Denote by A(D) the skew- 
symmetric matrix with the rows and the columns indexed by V, where a uv — a w ^ u,v ^ in case (u, v) is an arc 
of D, a uv = — a w ( u ' v } in case (v, u) is an arc of D, and a u , v — otherwise. 

Definition 4.5. The Pfafjian is defined as 

Pf G {D,a) = s ^s*{P)a lljl ■ ■•a in j n) 
p 

where P — {{iiji}, • • • , {injn}} is a partition of the set {1, . . . , 2n} into pairs, ik < jk for k = 1, . . . , n, and 
s*(P) equals the sign of the permutation iijx . . . i n j n of 12 . . . (2n). 

Each nonzero term of the expansion of the Pfaffian equals a w ^ or — where P is a perfect matching 

of G. If s(D, P) denote the sign of the term a w ^ in the expansion, we may write 

P/ G (fl,a)=^ S (APK< p ) 
p 
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4.2. Cubic lattices as generalized g-graphs. In this subsection we will describe how to realize Q{mx , k) 
as a subgraph of a generalized g-graph. Let us denote by Q L the larger lattice Q(mi, m,2, 2k + 1). 

How to draw Q L on the plane. First draw the paths V xy along a cycle in the linear order as in Def. 11.11 

Next, draw the horizontal edges inside this cycle, and the width edges outside of this cycle as depicted in 
Fig- El below where Q L = Q(3, 3, 5) is properly drawn. The figure also shows how Q = Q(3, 3, 2) is embeded 
in Q L = (3,3,5). 

We keep the following rule: the interiors of the curves representing h xyz and h( x+1 \ yz (w xyz and w x / y+ i\ z 
respectively) intersect if and only if z is even. We denote by C(e) the curve representing edge e. 




Figure 3. Drawing of a cube Q L = Q(3, 3, 3) in a plane. 

Now we modify Q L into a generalized g-graph Q' . We introduce new vertices to some edge-crossings 
of Q L and delete and subdivide some edges. In this way we obtain a generalised g-graph Q' with an even 
subdivision of Q(mi, JTI2, k) as its subgraph. The construction is analogous and independent for each plane, 
and so we describe it only for a plane W x . 

The construction is described by Fig. 01 for the edges between V x i v -x\,V X y and between V xy ,V x ( y+ i), for 
x odd and y < mi — 1 even. 




Figure 4. Construction of Q' for part of the plane W x (x odd), V x < y -i), V X y, V x ( y +i) (y even). 

(1) For each y even let Aux\ = {w xyz ; z odd }. For each edge e of Auxi introduce a new vertex to each 
intersection of C(e) with the curves representing the edges of W x / y -i) U W, where W = W x (y+i) in 
case y < m — 1 and W = otherwise. By this operation, each e 6 Aux\ is replaced by a path. Call 
each edge of this path auxiliary. 
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(2) For each y even let Auxi — {w x ( y — 1)1, w x (y— i)n} U A t where A — w x(y+i)n} hi case 
y < m — 1 and ^4 = otherwise. For each edge e of Aux2 introduce a new vertex to each intersection 
of C(e) with the curves representing the edges of W xy . Hence each e £ Aux2 is replaced by a 
path. Call each edge of this path auxiliary. For each y even the edges v xy i, v xy / n _x) an d also 
v x(y+i)ii v x(y+i)(n-i) wm also be called auxiliary. In Figure 01 the auxiliary edges are represented 
by dashed lines. 

(3) The edges w xyz , y even and z even will be called relevant for Q. If y < m — 1 then the relevant edges 
are subdivided by two vertices (added in 2.) into three edges of Q' . The middle one will be called 
special and the other two long. If y = m — 1 then the relevant edge w xyz is subdivided by one vertex 
into two edges of Q' . The one incident to V xm will be called special and the other one long. If e is a 
relevant edge of Q, then we choose a initial long edge / and we let w(e) = w(f). We let the weight 
of the special edge and of the remaining long edge be equal to 0. 

(4) The edges of W x ^ y _i- ) U W also got subdivided by new vertices introduced in step 1 and step 2. 

(5) We delete all edges of the paths obtained from w x f y _u z and w x f y +i) z , 1 < z < n odd, as well as we 
delete all vertexes in intersections of those edges with auxiliary edges. In Figure^ the deleted edges 
are represented by dotted lines. 

(6) Each edge e £ {w x ( y _i) z , w x ( y+ i) z ; z even }, is subdivided by new vertices introduced in step 1 into 
a path. We let the weights assigned to the edges of the path equal except of one initial edge whose 
weight is let equal w(e). The edge e of this path such that the interior of C(e) does not intersect 
interior of any curve representing a long edge will also be special. The others are called short. 

(7) All vertical edges which are not auxiliary (see FigQJ will be called special. In Figure 01 the special 
edges are represented by normal lines. 

(8) Let Aux denote the set of all auxiliary edges. Then deletion of Aux results in a subdivision of 
Qm 1 m 2 k- We subdivide some special edges (vertical and border width ones) so that the graph 
Q = Q 1 — Aux is an even subdivision of Q mi m 2 k- All these new edges will be special, and we set 
their weights equal 0. 

This finishes the construction. In Figure the long and short edges are represented by fat lines. 

Q' is a generalised g-graph. Its planar part Q p is set of all the auxiliary and special edges. Other 
edges (i.e. the short and long) are drawn on a face of Q p and they may be drawn onto a pair of bridges 
above this face. One bridge contains one long edge, and the other bridge contains all the short edges. 

Now let us set appropriate weights to all the auxiliary edges. In order to do that first subdivide and add 
some auxiliary edges and vertices to Q p in such a way that there exists a matching Ma U x of V(Q') — V(Q) 
consisting of auxiliary edges, and a perfect matching M of Q p . We set weights of edges from Maux as 
10(e) = 0. Weights of other auxiliary edges w(e) — —00. With such properties of Q' we have an important 
result: V(Q',a) = V(Q(m\, raj, k), a). 

Now we introduce a basic orientation D p of Q p in order to prepare our stage to use the Theorem 0] for 
the generalised g-graph Q' . Orientation D p has the following properties: 

1. D p on special edges is in agreement with the natural ordering (Def. II. 

2. All the signs s{D p ,P) in Pfaffian Pf QP {D p 1 a) (Def. 02) are positive. 

3. The orientation of edges on a bridge has positive sign if and only if it is in agreement with the natural 
ordering fDcf. II. if . 

The construction of such a D p is possible according to Kasteleyn [K] . 

5. The proofs 

5.1. Proof of Theorem ^ We use Theorem 0] for the graph Q'. Then we may write 

49 
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2=1 
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where the first equality is the result of previous section. The second equality holds because of the Theorem 
0| The third equality holds because of the setting of the weights of edges in Q 1 in the previous section. And 
finally, orientation D* being induced on Q by Di on Q, the fourth equality follows from the three facts: 
setting of the weights of edges, Q is even subdivision of Q, it is possible to permute vertexes covered by 
Maux so that the sign is correct. 

In Section |21 we have defined for a bipartite graph G a matrix Z{D) depending on a orientation D. It is 
easy to check that for Q holds Pf Q (D*,a) = det(Z(D*)). 

To proceed in proof of Theorem we have to establish the dependence of c(r on the orientation D* 
of Q and specify the orientations we sum over in the language of Q not Q' . 

Each relevant orientation Di of Q' is determined by the fixed basic orientation D p of Q p , and by a pair 
of signs for each pair of bridges. Each pair of bridges is associated with one long edge and set of short edges 
of Q'. Hence these signs may be given by specifying (dp. (e), <P D . (e)) 6 {H — } 2 , for each long edge e, where 
d l D . (e) denotes the sign of the bridge containing e, and d 2 D . (e) denotes the sign of the other bridge containing 
the set of short edges. 

The relevant edges w x i m2 -i) z {Q{mim2k)) and h( mi -\) yz {Q{mim2k)) are associated with only one long 
edge of Q'. If e is such relevant edge of Q, we will call it border edge, and we denote by e\ the corresponding- 
long edge. We let do^e) = (d]-, (e{),d 2 Di (ei),+,+). Each relevant non-border edge e of Q has two long 
edges e%, e2 associated with it. We let djj i (e) = [d D , (ei), d 2 D . (ei), dr D . (e2), d 2 D . (e2))). A relevant vector is any 
element r of [{+, — } 4 ] K such that r(e)z — r(e)± = + for each relevant border edge e of Q. 

From the definition of c(r(Di)) and from properties of the generalized g-graph Q' of genus g follows that 
c(r(Di)) = 2~ 9 (— 1)™, where n is number of faces of Q p above them edges on both bridges are oriented 
in negative sense, i.e. against fixed ordering. There is a natural bijection between relevant orientations 
of Q' and relevant vectors. If r is a relevant vector, then let Di(r) denote the corresponding relevant 
orientation of Q' and let sgn(r) of a relevant vector r be calculated as follows: sgn(r) = 2 9 c(r(Di)) = 

»=0,1; r(e) 2 » + i=r(e) 2j+ 2=-l}| i 

There is C r = km\ (rn% — l)/2 + kui2{m\ — l)/2 relevant edges and C& = k(m\ + ma) border edges. Genus 
of Q' is therefore g = 2C r — C&. There are 4 2C ^- C >' relevant vectors. 

Let us consider a set T> of relevant orientations Di such that there is at least one relevant non-border edge 
e with r(e)a ^ f(e)4, i-e. orientation of bridges with short edges incident to e are different. 

Proposition 5.1. Sum of contribution of all orientations from T> to V(Q, a) is zero. 

Proof: For given orientation D e V define a set Ajj = {e relevant, r(e) 2 7^ r(e)^}. Now let fix a set o 
relevant edges Af and consider all orientations Df such that Ajy. = Af. Fix also one relevant edge e/ G Af. 

Remark: c(r(Di)) + ciriD^)) = if r{D\) and r{D2) differs only in one relevant edge e in such a way 
that r Dl (e) 2 = r D . 2 (e) 2 ^ r Dl (e) 4 = r D2 (e) 4: and r Dl (e)i = -r D . 2 (e)i, r Dl (e) 3 = ~r D2 (e) 3 . It is also that 
induced orientations DI — D\. 

Using this remark for e = ej we have for every Af. 

]T c(r(D l ))Pf Q (D*,a) = 0. 

The claim in Proposition 15 . II follows directly. □ 

Let us call all Di (I T> useful orientations and from now we sum only over such orientations. There is 
2 3Cr ~ 2Cb of useful orientations. Similarly we define useful relevant vectors. 

Proposition 5.2. Each orientation D of Q, which agrees with the fixed order on vertical edges, corresponds 
to 2 Cr ~ Cb of useful orientations Di of Q' and all those Di have the same sign c(r(Di)). 

Proof: Let r be a useful vector. Then D(r) determines uniquely r(e)2 and r(e)i for each relevant edge e and 
also r(/)i for each relevant border edge /. Hence D(r) determines uniquely r(/) for each relevant border 
edge /. Moreover D(r) determines uniquely the product r(e)i x r(e) 3 for each relevant non-border edge e. 
Since there are C r — Cb relevant non-border edges, there is 2 Cr ~ Cb orientations Di corresponding to it. 

Let r, s be useful and both lead to the same D of Q. Then r(e)2 = s(e)2 = s(e)4 = r(e)i for each relevant 
non-border edge e and r(e)2 = s(e)2 and s(e)i = r(e)x for each relevant border edge. This implies that 
sgn(r) = sgn(s). □ 
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Concluding from what we did till now we can write the dimer partition function as 



V(Q, a) = 2~ c " sgn(D)det(Z(D)), 

D 

where the sum is over all orientations of Q with all vertical edges positive. At this point the only thing 
missing in proof of Theorem^is to concretize sgn(D). 

Proposition 5.3. sgn(D) = (-l)l{( e > e ')efl(<2); 8 ign(r», e )= s ign(£., e ')=-l}| ) where the set j edges #(q) ls defined 
inTUk 

Proof: Let us remind that sgn(£>) = (—1)" where n is number of faces with both bridges above them 
negative. Each non-border relevant edge e corresponds to two faces with bridges. The contribution of 
the two corresponding faces to sign could be negative only if orientation of e is negative (r(e)i ^ r(e)a). 
And moreover if orientation of corresponding shot edges is negative (r(e)2 = f{e)i = — !)• It follows from 
the construction that this happen only if product of orientation of all e' (as defined in 11.21 ) is negative. 
Considering all the relevant edges we get the claim of Proposition |^| □ 

5.2. Proof of Theorem HJ 

Proposition 5.4. The average ofdet(Z(D)) over all orientations of Q with all vertical edges fixed positive 
equals a w ^ M \ where M is perfect matching of Q consisting only from vertical edges. 

Proof: By the linearity of expectation the contribution of other than vertical edges cancel out when we 
calculate the average of det(Z(D)). Since Q(mim2k) has exactly one perfect matching consisting of vertical 
edges only, Proposition 15 . 41 follows . □ 

There is 2 2C orientation D in^to be summed over. And 2 C ~ 1 (2 C + 1) of them have positive sign. Proof: 
Looking at Figure |21 we can see that as soon as there is at least one relevant (shaded) edge negative there is 
so many orientation with positive as with negative sign. Only if all the relevant edges are positive the sign 
is also positive. It means there is 2 C more positive orientation then negative. 

By Theorem ^ an d Proposition 15.41 we have then 

V(Q, a) = 2- c [-2 2C a w ^ + 2£{2 C -\2 C + 1))] = -2 c a w(M ^ + £{2 C + 1), 
where £ equals the average of det(Z(D)) over all D of Q with positive sign. And Theorem|21is proved. □ 

5.3. Proof of Theorem |3J We use a straightforward reformulation of a theorem of |L2j : 
Theorem 5. Let G be a generalized g- graph. Then 

4 s 

H(G,a)=^c(r(A))n( 1 + W "(H)' 

»=1 Di 

where Di are all the relevant orientations of G, and rotation of p and I1_d (I + W([pD are defined in the 
same way as described for Q in section\^ 

Proof of this theorem could be found in |L2j and for planar graphs in [S] . 

Let Di, i — 1, . . . , 4 s be the relevant orientations of generalised g-graph Q 1 . If we apply theorem [S] to Q', 
from the choice of weights we get 

43 

i=l D* 

The same analysis as in the previous section finishes the proof of Theorem [31 
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6. Conclusions 



We express the finite 3D Dimer partition function as a linear combination of determinants of oriented 
adjacency matrices in Theorem ^ and the finite 3D Ising partition function as a linear combination of 
products over aperiodic closed walks in Theorem [3] 

The methodology we use is embedding of cubic lattice on 2D surfaces of large genus. By this embedding 
cubic lattice is transformed to a generalized g-graph. Then we use theorems from |GLIIL2| . which generalize 
approach of [K|. 

We believe this approach may relate to recent developments in conformal field theory |MJ ICM1I ICM2| 
and topological string theory [DRV , and ultimately teach us something about 3D basic statistical physics 
problems. 

References 

[ORV] A. Okounkov, N. Reshetikhin, C. Vafa, Quantum Calabi-Yau and Classical Crystals, hep-th/0309208 
[AMV] Alvarez-Gaume L., Moore G., Vafa C, Theta- functions, Modular invariance, and Strings, Commun. Math. Phys. 106 
(1): 1-40 1986. 

[CM1] Costa-Santos R., McCoy M.B., Dimers and the Critical Ising Model on lattices of genus > 1, Nuclear Physics B 623 
(3): 439-473 FEB 25 2002. 

[CM2] R. Costa-Santos, B.M. McCoy, Finite size corrections for the Ising model on higher genus triangular lattices, Journal 

of Statistical Physics 112 (5-6): 889-920 SEP 2003. 
[M] Mercat C, Discrete Riemann surfaces and the Ising model, Commun. Math. Phys. 218 (1): 177-216 APR 2001. 

[S] S. Sherman, Combinatorial Aspects of the Ising Model of Ferromagnetism I, J. Math. Phys. 1, 1960. 

[D] Distlcr J., A note on the 3D Ising model as a string theory, Nucl. Phys.B 388 (3): 648-670 DEC 28 1992. 

[LI] Loebl M., On the Dimer Problem and the Ising Problem in Finite 3-dimensional Lattices, Electr. J. Comb. R30 9 

(1): 2002. 

[L2] Loebl M., A Discrete Non-Pfaffian Approach to the Ising Problem, DIMACS, Series in Discrete Math and Theor. 

Comp. Science, 63 2004. 

[GL] A. Galluccio, M. Loebl, A Theory of Pfaffian Orientations I, II, Electr. J. Comb. 6 (1): 1999. 
[K] P. W. Kasteleyn, The statistics of dimers on a lattice, Physics 27: 1209-1225, 1961. 

Dept. of Applied Mathematics and, Institute of Theoretical Computer Science (ITI). Charles University, 
Malostranske n. 25, 118 00 Praha 1, Czech Republic. 
E-mail address: loeblSkam.mff.cuni.cz 

Institute of Physics, Academy of Science of the Czech Republic, Na Slovance 2, 182 21 Praha 8, Czech Republic. 
E-mail address: zdeborl@fzu.cz 



10 



